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INTRODUCTION ; 


Estimation of fluxes from satellite radiance measurements 
is known to be sensitive to the models used in the conversion and 
to the proper identification of the scene (Ruff et al., 1968; 
Raschke et al., 1973) particularly cloud. The models used in the 
radiance to flux conversion are the bidirectional and directional 
models for clear and cloudy (with different levels of cloudiness) 
scenes. The cloud identification for this purpose is derived 
from the narrow band 11.5pm channel on the Temperature Humidity 
Infrared Radiometer (THIR) on Nimbus-7 (Chen et al., 1980; Hwang, 
1982). The scene identification in the data processing stream 
for Nimbus-7 Earth Radiation Budget (ERB) instrument is based on 
bispectra 1 thresholds (Jacobowitz et al., 1984) on the broad band 
radiance measurements in the short wave (0.2 to 4.8 p m) and the 
long wave (4.5 pm to 50pm). Scene determination using thres- 
holding methods is shown to be deficient and some alternatives 
"ere suggested in Vemury et al . (1984). Data from the Earth 

Radiation Budget Experiment (ERBE) will be processed using a 
statistical procedure known as the Maximum Likelihood Estimation 

(MLE), for detection of cloud amount using the SW and LW radiance 
measurements . 

This scene information derived from NFOV channels is also 
used xn the inversion process for the medium field of view and 
wide field of view (MFOV and WFOV) channels on the ERBE instru- 
ment. Correct scene identification from the narrow field of view 
(NFOV) channels, is thus of great importance. 
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This report is based upon a study to test the MLE proce- 
dure and evalute its performance on actually observed radiance 
data from the narrow field of view (NFOV) channels on the Nimbus- 
7 ERB instrument. The Earth radiation budget parameters obtained 
with the MLE procedure will be compared with those obtained from 
the Sorting into Angular Bins (SAB) method. The present study 
also serves as a test of performance of the scene identification 
algorithm to be used in the ERBE data stream. 

2. Description of SAB and MLE Methods : 

Since scene selection is an identified problem with most 
inversion methods, the results obtained with the MLE method will 
be validated by using a data set which is independent of scene 
identification and corresponding scene dependent bidirectional 
models. Sorting into Angular Bins (SAB) method provides such a 
validation data set. (For details see Arking et al. f 1984). This 
method essentially involves the binning of the SW and LW radiance 
measurements into different solid angle bins (85 in all) for each 
equal surface area region on the globe referred to as a target 
area (TA). One month of data is sorted into these bins to provide 
adequate sample size in each of the 85 bins for each TA . The 
radiance means for each bin for the month are computed and for 
each TA, the 85 bin radiance means are integrated with the appro- 
priate solid angle weight for the bin to obtain the TA flux. 
This process involves no scene identification and is really 
independent of the underlying scene, viz., cloud or no cloud and 
also of the associated angular models. 
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A directional model is necessary to obtain a mean daily 
albedo from the instantaneous albedo for the TA . In the present 
study, mean daily albedos are not calculated in order to reduce 
the complexity of the problem. All the albedos referred to in 
the study are instantaneous albedos. 

The input data set for this method is the Sub-target Radi- 
ance Tape (STRT) which contains 9 days of data compacted into one 
tape. Some necessary definitions are given in Frame 1 and Frame 
2 and a very brief description of the SAB procedure is shown in 
Frame 3. On the same frame, a few comments are made on the data 
quality checks and processing methods. 

The binning procedure will fail if the sample sizes for 

the bins are not adequate so that one is required to fill a 

number of bins with interpolation for the TAs . The number of TAs 

j d due to lack of adequate sample sizes is found to be 

rather small in all cases, viz., reflected radiances, emitted day 

time and emitted night time radiances. Figures 1 and 2 show this 

result schematically on a latitude band basis. The abscissa is 

the zone number, 1 for the South Pole and 40 for the North Pole; 

each latitude band is 4.5' wide. At each latitude, the number of 

TAs rejected is shown by the dashed line while the solid line 

indicates the total number of TAs in that band. It is obvious 

that a very small number of TAs are being rejected with the 
binning process. 

The SAB algorithm consists of software at three levels. 

At Level I, the radiance observations are grouped into TAs and 
for each TA , into angular bins and radiances and sample sizes are 
accumulated for nil the days in the month. The Level II algo- 
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rithm uses the output from Level I stored on tape and after 
quality control checks, fills any empty bins and for each TA , if 
accepted, integrates the radiance means over the solid angles to 
obtain the flux. Output products for Level II are the reflected 
and emitted fluxes (day and night) for each TA with a few subsid- 
iary pieces of information. Level III software computes the 
zonally , and globally averaged earth radiation budget parameters. 

Maximum Likelihood Estimation : 

The input data set for this method is the Master Archival 
Tape (MAT) which consists of corrected earth located radiance 
values. There is one MAT tape for day. The other input is the 
table of anisotropic factors (reflectance and emission) for each 
of 12 earth scenes identified for use by ERBE algorithm. A brief 
description of the MLE procedure is given in Frames 4 and 5. No 
further details about the MLE method will be discussed since the 
ERBE team is quite familiar with most of the details. 

3 • Operational Steps 

The following sequential steps were performed to retain 

the integrity of the original MATRIX algorithm, while the modifi- 
cations were made. 

(i) The MATRIX algorithms were modified to be able to 
process the NFOV data only on a daily basis. The 
results of a few days of data were compared with 
MATRIX output. 

(ii) The code was further modified to compute the 
instantaneous albedoes only. No directional 
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models were used in the present study. 

(iii)(a) The MLE scene identification algorithm was set 
up for latitude, longitude directions different 
from those accepted by Nimbus-7 ERB MATRIX 
algorithms. Also, the spatial resolution and 
global grid structure used by the MLE system are 
different. Therefore the Scene ID was performed 
on the ERBE grid structure, but the derived fluxes 
were stored for the MATRIX TA grid. This is 
necessary to make the comparison with SAB 
final output product. 

(b) The relative azimuth used by ERB is the complement 
of the relative azimuth as used by MLE. 

(c) Underlying scene is the climatology scene used by 
ERBE system on the 2.5° X 2.5° grid instead of the 
target area system used by Nimbus-7. 

(d) The bidirectional models for reflectance and the 
emission models are the ERBE models. 
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(e) Data handling and rejection of non-valid observa- 
tions are different and dealt with at different 
locations in the code in the ERBE system. Some 
modifications were made to the original Langley/ 
NASA supplied code to handle these appropriately. 

(iv) Several test runs on the data revealed a few 

computer-systems related problems. These were 
identified and corrected. 

(v) Observations were truncated at different satellite 
zenith angles from 90° to 40° and the effect on 
the daily averaged instantaneous albedo, LW (day, 
night and total ) fluxes was investigated. 

Steps (vi) and (vii) were performed, using the emission 
models and bidirectional reflectance (ERBE) models from Nimbus-7 
and directional models developed by ERBE working groups (ERBE 
unpublished document, 'Angular Radiation Models for ERBE': 
October 30, 1984). 

(vi) Process each of 23 days in June 1979 through the 
algorithms and produce the following data sets: 

(a) Average instantaneous albedo and LW flux 
(day, night, total) for each of the 2070 
target areas on the globe. 

(b) Different levels of cloudiness percentages 
also for each TA . 
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(c) 

(d) 

(e) 
( vii) 

( viii ) 

(ix) 

(x) 


All the individual converted flux and 
albedo values for chosen (20) TAs. 

All information regarding viewing and 
incident angles and scene selected, SW 

and LW radiances and the computed fluxes 
for 20 TAs. 

Statistical parameters which would provide 
information on the model confidence levels 
for all 12 surfaces. 

From (vi) above, monthly averages of the 
parameters in (a) and (b) under (vi) for 
each of the TAs were computed. 

SAB: From Level I data set of the SAB 

method. Level II algorithms were applied 
to create the TA averages of the instan- 
taneous albedo and the LW fluxes. From the 
TA averages, the zonal and global mean values 
were computed. No Level III algorithms were 
used in the present study. 

The TA zonal and global means of the param- 
eters from SAB, and from the use of Nimbus-7 
models were compared. 

Cloudiness percentages were converted from 
TA to zonal averages for each cloudiness 
condition from clear to completely cloudly 
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(xi) 


(xii) 


4 . Results : 

(Results of the validation study were presented at the 
ERBE Team Meeting at NASA/Langley Research Center on November 1, 
1984.) 

In this section, we discuss several of the results which 
may be classified broadly as: 

(i) the final monthly mean ERB parameters 

(ii) the monthly mean cloudiness - dependence on 
latitude band 

( i ^ ) the effect of eliminating all radiance observa- 
tions beyond a satellite zenith angle 

(iv) Scene Identification Reliability Index. 

Several of these items have come up as offshoots of the original 
purpose of the study. Some of the available data and results 

pertaining to the 20 TAs mentioned in Section 3 (vi) C and 


for each day. The percentages were plotted. 

Monthly mean cloudiness percentages were com- 
puted in both cases from the daily means. 
Zonally averaged plots of the monthly mean 
cloudinesses are produced. 

Similar plots were also produced for each of 
the studies with truncation of observations 
at different satellite zenith angles. 
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3 (vi) d have not been analyzed due to limited funding under the 
present contract. 

(i) Monthly Mean Earth Radiation Budget Parameters : 

These parameters are different from the conventional 
parameters used in most radiation budget calculations. For ex- 
ample, the albedo, computed and presented here is the instantane- 
ous albedo corrected only for anisotropy of the reflecting sur- 
face and a diurnal correction is not applied. The LW fluxes are 
corrected for the limb brightening or darkening as the case may 
be. The LW flux is assumed to be constant through the night or 
the day, and no models are applied for the time-dependent varia- 
tions. The monthly means are obtained from the daily means as 
simple average over the number of days used in processing. In the 
present case, data for 1979 are used and there are 23 days of 
observed data from the Nimbus-7 ERB instrument. The computed 
monthly mean values are shown in Table 1 as global averages. Net 
radiation is not obtained because the albedo presented here is 
not the daily mean albedo. For purposes of comparison with the 

validation data set, similar values obtained with SAB are also 
shown . 

We notice from Table 1 that the global mean radiation 
budget parameters of interest to us, computed with the MLE method 
show good agreement with those obtained through the binning 
procedure using SAB method. Keeping all observations upto 90° in 
satellite zenith, the instantaneous albedo is 28.1 with the MLE 
method while the SAB value is 27.4. Using the bispectral thres- 
hold method applied in Nimbus-7 data processing, the monthly mean 
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diurnally corrected global albedo is 33.0 with a corresponding 
value of 30.1 from SAB method (Arking et al., 1984). Of the three 
ERB parameters, albedo is the one which showed worst agreement 
between SAB and MATRIX NFOV product. This difference of nearly 
2.9 in albedo is now reduced to 0.7 with the MLE method. The 
improvement gets better when observations are truncated at 75° 
in satellite zenith angle. The LW fluxes (day time, night time 
and total) also improved with this truncation. The day time 
longwave flux difference with SAB is 2.7 w/m 2 with the SAB value 
at 243.7 w/m 2 and MLE value of 241.0 w/m 2 . Night time and total 
flux values show even better agreement with SAB values. The 
effect of truncation at smaller angles is not appropriate and the 
reason will be discussed in a later section. 

Differences in the zonal averages using the SAB and MLE 
with cutoff at 90° are shown in Fig. 3a for the albedo, and Figs. 
4a, 5a, and 6a for the day time LW flux, night-time and total LW 
fluxes for the monthly averages respectively. On a zonal basis, 
the instantaneous albedoes with MLE are slightly larger by an 
amount smaller than 0.015 and for some latitudes in very good 
agreement with the SAB zonal averages. The day time LW fluxes 
from SAB are about 4 to 6 w/m 2 larger on an average, while the 
night-time values are different by about 2 watts in the same di- 
rection. The zonally averaged total LW flux is also larger 
(Fig. 6a) by about 3 w/m 2 . 

The TA fluxes with MLE are compared with those using SAB 
in Figures 3b, 4b, 5b and 6b. The SAB values are shown on the 
abscissa and the MLE values on the ordinate. The statistical 
parameters related to the linear regression are shown in Table 2. 
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Concentrating for a moment on the instantaneous albedoes, we no- 
tice that in the large albedo region, there is a good agreement 
between the SAB and MLE albedoes. In the low albedo region 

(i.e., between 0.1 and 0.3), MLE tends to overestimate the al- 
bedo by as much as 0.1 in some cases. In the intermediate albedo 

range, the opposite may be the case, with inaccuracies of the 
order of 0.05 or less. 


We could conclude that in the completely cloudy and clear 
snow cases (high albedo), the angular models appear to behave 
very well while in situations of partial cloudiness as also in 
some clear cases, the errors vary between 0.05 and 0.1. These 
conclusions, however, need to be confirmed with individual albedo 
computations for some chosen TAs, before a final definitive 
conclusion can be drawn. This comparison is in progress. 

On an overall basis, the albedoes with MLE are regressed 
against the SAB albedoes and the gradient, intercept and the co- 
efficient of regression assuming a linear relationship are shown 
on the figure. The standard error of regression is 0.015. Value 
of o e , where o q 2 is the sum of the squares of the deviations 
between MLE and SAB values, is also shown on the figure. 

Day time fluxes are shown in Figure 4b. The parameters of 
regression are also shown. The agreement at large fluxes is much 
better between the two methods, corresponding to warmer or non- 
cloudly land or desert areas. The standard error from the linear 
regression is about 4.5 w/m 2 . 

In the night time, similar results apply (Fig. 5b) except 
that the agreement is better in very low LW flux region. These 
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apparently correspond to observations from TAs which have night 
time all the 24 hours (or from the south polar regions and the 
low fluxes confirm that). From Figures 4b and 5b, we may observe 
that the polar long wave models are performing quite well. On a 
global basis, the night time LW fluxes show a smaller coefficient 
of regression compared to the day time or total fluxes. Similar 
comparisons for the total flux are shown in Fig. 6b. 

Effect of eliminating observations greater than 75° is 
shown in Figures 7a, 8a, 9a and 10a. The zonal averages are 
shown as the difference between the value with SAB and the albedo 
with MLE when they are truncated at 75°. The regression plots are 
shown in the corresponding Figures 7b, 8b, 9b and 10b. The effect 
of truncating observations at 75° is to decrease the zonal 
averages of albedo generally over the southern hemisphere, and 
the effect is much smaller over the northern hemisphere. The LW 
fluxes seem to be enhanced when the observations at large satel- 
lite zenith angles (>75°) are removed and this is true in all the 
LW flux computations (day, night and total). In all cases, the 
SAB, Nimbus-75 differences are smaller than with the Nimbus-90 
case . 

Mercator projections of the albedo and LW flux fields 
between 70° S and 70°N are shown in the next 12 figures (i.e., 
Fig. 11 through Fig. 22. Fig. 11 through 14 are the parameter 
fields with the SAB method. As the earlier results showed better 
agreement between SAB and MLE with cutoff of 75°, the MLE results 
with 75° cutoff case are shown in Figs. 15 — 18. The difference 
plots are shown in Figs. 19-22. 
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( i i ) Monthly Me_an_ C_lo_ud > P^ej"_c_eji_tj|_g_es. * 

The presence or absence of a cloud and in case, a cloud is 
present, the degree of cloudiness is determined by the scene al- 
gorithm using the MLE probability density function. The percent 
cloud under each category of cloudiness is determined from the 
total number of observations of the TA and the number of cloud 
identifications in each category. The TA cloudiness categories 
are then composited to obtain the zonally averaged cloudiness 
percentages under each category for each day and the zonal mon- 
thly mean cloudinesses are computed. 

Figure 23 shows the monthly mean zonal average cloudiness 
on the basis of the MLE method. The scene detection algorithm 
uses LW radiances only in the south polar regions due to the non- 
availability of SW radiances. All the SW and LW radiance 
combinations at the North Pole indicate a high degree of cloudi- 
ness near 100%. The position of the inter-tropical convergence 
zone seems to be very well defined slightly north of the equator 
and a region of high cloudiness at the South Pole. In addition, 
a regional peak of cloudiness in the southern hemisphere mid- 
latitude and a small peak and a plateau of cloudiness in the nor- 
thern hemisphere mid-latitude can be observed. These may be 

typical of the convergence zones of the ascending portion of the 
Hadley cell. 

Figure 24 shows the percentage clear for each latitude 
band in a similar way. In contrast with the cloudiness figure, 
not even one latitude band of the globe is completely clear. A 
maximum of near 55% of the zone classified as clear appears 
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around 40°N in June, with two peaks almost of equal percent clear 
in the equatorial north and south (around 20° N and 20° S respec- 
tively). The region of ITCZ clearly marked in Fig. 23 is almost 
0% clear from Fig. 24 indicating that the remaining 50% of the 
latitude band is either partly or mostly cloudy. 

Figures 25 and 26 are for the partly cloudy and mostly 
cloudy scene identifications and thus, complete the monthly mean 
picture. The winter hemisphere has greater degree of cloudiness 
in both these figures and in no case, do they exceed 40% of the 
latitude band. Combining Figures 23, 25 and 26, we would note 
that the degree of cloudiness is much larger in the winter hemis- 
phere. Such a conclusion can also be drawn from looking at clear 
percent plot of the southern vs northern hemisphere in Figure 24. 

(iii) Effect of Satellite Zenith Angle Cutoff : 

(a) Budget Parameters : 

Many studies analyzing the data from satellites have 
noted the ambiguity of scene selection at large satellite zenith 
(Raschke et al., 1973; Vemury et al., 1984). The effect of trun- 
cating observations at different satellite zenith angles on the 
albedo and LW flux have thus been investigated as a function of 
the satellite zenith angle. Table 3 provides a summary of the 
global mean values (with standard deviations) of the present set 
of budget parameters. They are also shown in Figures 27 and 28. 
The effect was studied for a day (day 152 of 1979) which happens 
to be the second day of a three day on and one day off cycle and 
thus has the maximum number of observations. As expected from 
earlier studies, the global albedo drops as the satellite zenith 
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threshold is decreased. Also, to be noted is the difference in 
the albedo gradient for angles 90° to 70° and for 60° to 40°. Be- 
tween 70 and 60°, there is a sudden drop in the albedo. This is 
the region, where the number of observations by the scanner 
double at 58°, discussed in Vemury et al. (1984). 

It is also important to note that the effect of truncation 
of observations also has the effect of eliminating some TAs in 
contributing to the global mean. In the case of instantaneous 
albedo, for example, during one day 1849 TAs contributed to the 
global mean. With a total of 2070 equal area regions, nearly 90% 
of the TAs participated in providing the global mean albedo. The 
TA sampling drops rather drastically at 40° cutoff, with only 60% 
of the globe by area contributing to the global mean. The situa- 
tion becomes more severe for night time LW flux, where at 40 cut- 
off, 50% of the TAs provide a contribution to the global mean 
computation. In view of this reduced sampling with decreasing 
threshold value, it is necessary to choose a reasonable mean, 
where there is a compromise between the effect of incorrect scene 
selection and a poor sampling strategy of the global areas. The 
table also provides values of the LW fluxes and corresponding 
sample sizes for the different thresholds. As is apparent, lower 
than 70°, the number of TAs providing the global mean drops dras- 
tically due to the larger gradient. Thus removing too many ob- 
servations (or choosing a zenith angle threshold smaller than 
70°) is not a good sampling strategy. 

I-or comparison, the values with the Nimbus-7 ERB threshol- 
ding method are shown as dashed line in Fig. 27. The scale is on 
the left on the inside. These are the mean daily albedo values 
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(Vemury et al., 1984) with bispectral thresholds, with cutoff at 
different angles, for June 22, 1979. It is apparent that the 

thresholding method for scene selection leads to a much larger 
decrease in albedo at almost all cutoff angles. While a perfect 
scene selection and correct angular models should not show any 
slope with satellite zenith, the result with MLE indicates that 
there is considerable improvement over the bispectral threshol- 
ding method. 


(v) Scene Identification Reliability Index: 

A measure of performance of MLE as a scene selection 
scheme may be obtained by defining a parameter called the Scene 
Identification Reliability Index. Once a scene type is chosen 
based on the short wave and long wave radiance measurements, this 
parameter indicates to within how many standard deviations, the 
measured and model radiance values agree. The necessary 
definitions are shown in Frame 6. If both SW and LW radiances 
are available, scene selection makes use of both these. If only 
the LW or the SW is available, a scene is still selected. The 
percentage of observations which are within one standard devia- 
tion and two standard deviations under each of the categories are 
shown in Table 4. Nimbus-7 bidirectional and ERBE directional 
models (ERBE unpublished document, 'Angular Distribution Models 
for ERBE', October 30, 1984) are used in this case and the re- 
sults are shown for day 152 and day 153 in 1979. These percen- 
tages provide information on the level of performance of the mo- 
dels. For example, during the day-time for clear (0 to 5% cloud) 
ocean, there are 20,651 LW observations, of which 95.3% (under 
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NCLDAY in Table 4) are within two standard deviations of the 

model mean value. Of the 14,645 SW observations on the same day, 

almost all (99.99%) are within two standard deviations ( NCSDAY in 

Table 4). But only 94.6% of the observations with both LW and SW 

measurements, are within two standard deviations for both the SW 

and LW. During the night time, LW only data are used and for 

ocean, 98.1% of the observations are within two standard devia- 
tions. 

Behavior of clear snow seems to be opposite to that for 
ocean. During the day, 76.9% of the SW observations are within 
two standard deviations while nearly 96.9% of the day-time LW ob- 
servations are within two standard deviations. Almost, all models 
seem to perform well at night. When both SW and LW observations 
are used, snow models behave rather poorly during the day-time 
with only 38.5% of the observations (SW and LW together) within 
two standard deviations. Day-time desert observations seem to be 
in quite good agreement with model expectations for the SW (98.0% 
within two standard deviations), while for the long wave, the 
percentage drops to 68.7%. 

These different results for the 12 identified surfaces are 
shown in Fig. 29 for the case of the index less than or equal to 
2. For all the mostly cloudy scenes and for the overcast case, 
observations seem to be in very good agreement with the 
models at the 2 0 level, while the snow model seems to be the 
worst, followed by the desert model. Similar histogram of the 
model performance index at the 1 0 level is presented in Fig. 30. 

The "overall model" performance for the day is obtained 
from the bottom line of Table 4 for day 152. When both SW and LW 


17 



radiances are used, 90.3% of the 163,670 observations are within 
20 of the respective chosen models. During the night time, 
96.13% of the 229,946 observations are within 20 . 

Tho model reliability parameters thus indicate that the 
scene identification method is picking the scenes consistent with 
the angular model statistics. The better performance at night is 
not necessarily due to the models being more appropriate for 
night, but that the bidirectional model standard deviations for 
LW are much smaller than the SW standard deviations. 

5 . Conclusions 

1. Maximum Likelihood Estimation method provides 
albedoes and LW fluxes in very good agreement 
with the SAB method. 

2. MLE method has improved the scene selection 
compared to the bispectral threshold method. 

The albedo decrease (due to incorrect scene) 
is considerably reduced (Fig. 27). 

3. Adequate sampling strategy requires that 
satellite zenith angle cutoff should not be 
lower than 70° . 

4. There is still a remnant of the scene id. 
problem. Part of this may be due to the 
angular models as well. 
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5. 


As a test of the .software, the scene id. 
algorithm is performing well and there are 
no surprises. 

6. MLE procedure could pick out all scenes 
including partly cloud scenes. 


6. Future ERBE-Related Studies 

A. Comparison of ERB parameters on a regional 
or target area basis between SAB and MLE 
methods . 

B. Effect of GOES models (monthly). 

C. Effect of zenith angle cutoff with GOES 
models. 

D. Sensitivity of scene selection to LW-SW 
correlation coefficients. 

E. Are LW models adequate? 

F. Conduct experiments to determine the 
uncertainty in the flux measurements. 

G. Validate scene selection from other 
data sets. 

H. Continue similar studies for at least 
one more month. 
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Frames 


Definitions 


Inbtuntiinuouii solar irrudluncu measurement 

s„ W « 0 ) 2 

Solar zenith angle at the time of measurement 


Filtered radiance (measurement) 

Unfiltered radiance (after spectral correction). 
R= Reflected; E= Emitted 

Satellite zenith angle 


Relative azimuth 


Frame 1 



Bidirectional Reflectance: 


P(C') 


N (C*) 
S* cost' 


Directional Reflectance: 

2tt i 

r(C') » f [ p <C) sin0 d(sin0) d(j> 

■'o J 0 


Normalized Reflectance: 

r(Q 

r(C-O) " R(C) 


Instantaneous Albedo 

7TN r (C’) 
S’ C 08 C’ 


Angular Model Correction Factor: 


a (F’;e,<t») 


277 v/2 

{ | N r^' »®»^) sin ® d(sin0) d<J> 

0 0 


Instantaneous Corrected Albedo: 


TT N R (5* ) 
S' C085 1 


j A (c';0,<(») 


Frame 2 



SAB 


Bidirectional Reflectance 


P 

j,k,n 


TTN 


j,k,n 


S 0 (d/d o ) 2 cos 5 ' 


( 7 ) 


•3 “ 1» (observation number) 

k " 1. 85 (Bin 0 in TA n) 

n " 1. 2070 (TARGET AREA 0) 


f j" Fractional Field of View belonging to the TA n. 
Instantaneous Corrected Albedo for bin k 


m 



f P 

3 3> k »n 




( 8 ) 


(In this case f tis set equal to 1 ) 


“k ’ “k - cos0 s1 ” 9 « bln 

- ‘ 85 
Instantaneous Target Area Albedo « E A 

k»l k,n k 


(9) 


Frame 3 
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SAB 


1* tapes are used due to data compaction. 

2. Data quality is high due to severe quality control checks. 

3. Solar zenith angle information on the STR tapes is not 

precise. v 


4 . 


Fields 


of view which belong to 
by the proper fraction 


different TAs are apportioned 
of the radiance to each TA. 


5. Target Areas are rejected if 

a) the total sample size is less than 300. 

b) the center bin has less than 5 observations. 

c) any TA with more than 8 bins that have less, than 5 observations. 


6 . 


If a Target. Area is accepted, any empty bins are filled by 
interpolation from the near neighbors. 


7. The radiance means are computed and weighted by the solid 

angle of the bin and integrated over the 85 bins 
to obtain the TA flux. 


8 . 


The zonal and global averages are 
fluxes. 


calculated 


form the TA 


Frame 3 (Contd.) 
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maximum likelihood estimation 

REQUIRED SUBSIDIARY DATA SET 

From NIMBUS— 7 Scanner Database 

I. BIDIRECTIONAL MODELS (REFLECTANCE) 

(A) For each of 10 solar zenith angle ranges 

(B) For each of 12 scene types 

(C) For each of 49 angular bins 

(i) Obtain mean SW and LW ansotropic factors 

(ii) Obtain standard deviations 

(iii) Obtain correlation coefficients 

II. DIRECTIONAL MODELS (REFLECTANCE) 

(A) For each of 12 scene types 

(B) For 10 solar zenith intervals 

i) Albedo values 

ii) Normalized (to nadir) reflectances 

III. LONGWAVE MODELS (EMISSION) 

(A) For each of 10 Latitude bands 

(B) For each of 7 Satellite zenith intervals 

i) Amisotropic correction factors for the 
LW radiance to flux conversion 

ii) A priori Longwave mean values 

GEOGRAPHY DATASET 

A global dataset of surface classification on a 2 ♦ 5° x 2 • 5° 

into ocean, land, snow, desert and mixes land/ocean areas. 

Frame 4 


grid 
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ASSUMPTIONS FOR APPLYING MLE ; 

(A) That each set of SW, LW radiance measurements under 

each of the above cases is completely independent 
of other observations (statistical independence 
of -'the sample) 

(B) That the sample size is adequate and that the samples 

are distributed randomly to constitute a normal 
distribution 

(C) Complete, sample of SW, LW radiances belongs to a 

bivariate normal distribution. 


Frame 4 (Cont’d) 
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COMPUl AT ION : Probability density function under condition -CND is 


p.™ m“1 


CND 


"Qcnd 


2 " °CND< SU > °CND (LW) A - Pc 


P CND 


where 


P 

CND 

a 

P 

M 

L 


■ 

1 

[r-c] 

2 

i 

r - c) 

CND 2(1 - o 2 ) 

1 P CND J 

°CND (SW > 

► T \ 

°CND aW) 

M 

4 


• 


- 2 p 


CND 


is probability density 


a run ( SW ) a (LW) 
CND CND 


is one of four conditions of the viewed surface, 
viz., clear, partly cloudy, mostly cloudy or overcast 

is standard deviation from the model 

is correlation ceofficient 

Measured radiance values 

Apriori radiance values 


PROCEDURE : 


To identify scene, compute P for each of the 4 conditions for 
the viewed surface; pick the scene which has the maximum probability 
density function P. 


Frame 5 


28 



SCENE ID RELIABILITY CRITERIA 


LW DISP 


SW DISP 


(N 

- N 

) 

/ 



LW, 

MEAS 

LW, 

MODEL 

/ <r 

LW, 

MODEL 

(N 

- N 

) 

/ 



LW, 

MEAS 

SW, 

MODEL 

/ <r 

SW, 

MODEL 


NCLA ^ 2 
NCLDAY ^ 2 
NCLNIT ^ 2 
NCSDAY ^ 2 


LW . DISP ,<; 2 
LW . DISP ^ 2 
LW . DISP ^ 2 
SW . DISP 2 


AND SW. DISP ^ 2 
(DAYTIME LW DATA ONLY) 
(NIGHT LW DATA ONLY) 
(DAYTIME SW DATA ONLY) 


Frame 6 
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10. List of Tables 


Table 1 : 


Table 2 : 


Tab 1 e 3 : 


Ta b 1 e 4 : 


Monthly mean global averages of the radiation 

— n 86 m. ParameterS . COtnputed usin 8 i) SAB method 
ii) MLE method with Nimbus— 7 models with zenith 
angle cutoff at 90° and 75°. 

Notice that the agreement is better with SAB when 
Nimbus-7 models with cutoff at 75° are used. 


Regression relations for the SAB - Nimbus 
comparisons . 


Effect of zenith angle cutoff of observations 
on the instantaneous albedo, LW (Day, Night and 
Total) fluxes for Day 152 (June 1, 1979). 


The corresponding sample sizes are shown in 
Table 3B. 


Scene Reliability Index for each scene and 

cumulative value for the entire globe for 

Day 152 and Day 153. Nimbus— 7 models are used. 
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Table 1 


COMPARISON OF GLOBAL AVERAGES 
SAB AND MLE METHODS 


MONTHLY MEANS FOR JUNE 1979 



SAB 

MLE 



SAT. ZEN ^ 90° ** 

SAT. ZEN ^ 75° 

Inst . 

* 0.2738 

0.2813 

0.2804 

Albedo 




LW Flux 

243.7 

237.7 

241.0 

(Day) 




LW Flux 

232.8 

229.7 

232.2 

(Night) 




LW Flux 

235.6 

231.8 

234.7 

(Total) 





* With Diurnal Correction: NFOV (SAB) ALBEDO = 0.301 

NFOV (MATRIX) ALBEDO « 0.330 


** 


This case is referred to as the Nimbus-90 (N90) case through out the report to indieate that 
result s obtained with MLE with Nimbus-7 models with all observations up to 90° in satellite 
zenith angle. Under similar conditions when observations are limited to 75° in the zenith 
angle, the results are referred to as Nimbus-75 (N75) . 




A* B 

Inst. Albedo (N75) 0.017 0.960 

(N90) 0.018 0.968 

Day (LW) (N75) -0.032 0.995 

(N90) 0.715 0.978 

Night (LW) (N75) -0.889 1.002 

(N90) -1.196 0.992 

Total (LW) (N75) -2.044 1.005 

(N90) -1.290 0.989 


R 

a 2 

0.993 

0.2227 X 1CT 3 

0.992 

0.2402 X 10" 3 

0.998 

15.67 

0.997 

19.75 

0.997 

13.63 

0.996 

20.91 

0.997 

8.598 

0.995 

12.953 










Table 3. 


DSNAME=W5SKV.LIB.MLE 


DESCRIPTION NIMBUS 
J (90 DEG) 


DAILY MEAN 

ALBEDO 0.2812 


LONG WV.FLUXI 

(DAY) 1236.7007 


LONG WV.FLUXI 

(NIGHT) i 228. 126 A 


TOTAL 1230.5099 | 


(DLGLBTBL) 

_____ 000i0003 

. 00020007 

MLE 00030007 

00040007 

00050007 

DAILY 00060007 

GLOBAL AVERAGES 00070007 

JUNE 1ST 1979 00080007 

00090007 

I — | 1 00110007 

NIMBUS NIMBUS NIMBUS MTMnnc ...u.,,. 00120007 
(80 DEG) (70 DEG) (00 Si) ^"Si, gg||SSg? 

00150007 

00160007 

0.2802 0.2787 0.2727 0.2093 0.200 00U00O7 

00190007 

00200007 

00210007 

239.3793 241.5121 242.1,80 242.7481 243.5,5 ggflggg? 

00240007 

00250007 

00260007 

231.0549 232.5494 232.4121 232.4225 231.838 ggflggg? 

00290007 

00300007 

00310007 

233.3534 235.3143 235.7482 234.3254 234.447 ggllgggj 

! 100340007 

100350007 


Table 3(a). 


SAMPLE SIZES 






Table 4 


XXXXXXMXXKXXXXXKXXXXKKXKKXXXK 
* NIMBUS-7 x 

NXXHXMXXXXXXHXXXXXKXXXXXKXKXX 


xxxxxxxx day 152 xxxxxxxxxxx 


PRCNT 

.LE.X 


PRCNT | 
.IE. 2 | 


TOTAL PRCNT 
.LE.l 


PRCNT | 
.LE.2 | 


TOTAL PRCNT 
.LE.l 


PRCNT | 
LE.2 | 


TOTAL PRCNT 
.LE.l 


OCEAN 

LAND 

SNOW 

DESRT 

MIXED LAND/OCEAN 
PARTLY CLOUDY/OCEAN 
PARTLY CLOUDY/LAND 
PARTLY CLOUDY/ LAND# OCEAN MIX 
MOSTLY CLOUDY/OCEAN 
MOSTLY CLOUDY/LAND 
MOSTLY CLOUDY/LAND, OCEAN MIX 
COMPLETELY CLOUDY 


0.6682 
0.3792 
0.0000 
0 . 2490 
0.537* 
0.6659 
0 . *95* 
0.6110 
0.5022 
0 . **5* 
0.507* 
0.4821 


0.9*60 

0.8*36 

0.38*6 

0.65*2 

0.9222 

0.9806 

0.9*83 

0.9222 

0.9595 

0.9197 

0.9*99 

0.8820 


20651 

13557 

1051 

6669 

1838 

21509 

9593 

1059 

30982 

6**3 

1033 

107999 


0.71*6 

0.535* 

0.9163 

0.373* 

0.6028 

0.793* 

0.8985 

0.807* 

0.7069 

0.8006 

0.7202 

0.6797 


0.9529 

0.8655 

0.9686 

0.6866 

0.9255 

0.9858 

0.9756 

a.9*2* 

0.9710 

0.9781 

0.9632 

0.9202 


21278 

7665 

6110 

3631 

939 

2*503 

7132 

736 

37653 

8325 

762 

111212 


0.8067 

0.88*3 

0.6810 

0.7761 

0.8839 

0.9836 

0.9933 

0.9986 

0.8301 

0.9216 

0.9213 

0.6697 


0.9808 

0.9956 

0.9989 

0.9956 

0.9968 

1.0000 

1.0000 

1.0000 

0.9972 

0.9950 

1.0000 

0.9255 


0.9717 

0.7620 

0.1538 

0.6767 

0.9155 

0.9225 

0.5886 

0.7985 

0.7501 

0.6023 

0.7713 

0.6686 


xxxxxxxx DAY 153 xxxxxxxxxxx 



OCEAN 

LAND 

SNOW 

DESRT 

MIXED LAND/OCEAN 
PARTLY CLOUDY/OCEAN 
PARTLY CLOUDY/LAND 
PARTLY CL OUDY/L AND# OCEAN MIX 
MOSTLY CLOUDY/OCEAN 
MOSTLY CLOUDY/LAND 
MOSTLY CLOUDY/LAND# OCEAN MIX 
COMPLETELY CLOUDY 


0.6815 
0.3536 
0.0909 
0,2*97 
0.5668 
0.6*06 
0.53*1 
0.5316 
0 . *863 
0 . *885 
0 . *670 
0 . *828 


0.9*90 

0.8387 

0.8182 

0.69*9 

0.9515 

0.9798 

0.9572 

0.8921 

0.9598 

0.9*51 

0.9201 

0.8892 


TOTAL 

PRCNT 

.LE.l 

.L 

PRCNT | 
E.2 | 

TOTAL 

PRCNT 

.LE.l 

PRCNT 
.LE.2 | 

19627 

0.726* 

0 

.9569 

19967 

0.81*2 

0.9833 

13210 

0.5213 

0 

.85*8 

8900 

0.8522 

0.9967 

825 

0.9321 

0 

.9927 

*956 

0.5920 

1.0000 

6752 

0.3826 

0 

.7167 

3528 

0.8339 

0.9989 

1678 

0.651* 

0 

.9583 

1301 

0.81*8 

0.9992 

20950 

0.7633 

0 

.9850 

2*0*9 

0.9838 

1.0000 

87*1 

0.9119 

0 

.9815 

9039 

0.9933 

1.0000 

989 

0.7*62 

0 

.9171 

1013 

0.9990 

1.0000 

26979 

0.7101 

0 

.9723 

38*79 

0.8286 

0.9977 

6187 

0.8185 

0 

.9819 

8636 

0.9089 

0.995* 

9*0 

0.7021 

0 

.9500 

10*0 

0.8971 

1.0000 

99180 

0.6632 

0 

.9231 

109070 

0.6912 

0.9*65 

20605C 

0.68*0 

0 

.9328 

229978 

0.78*8 

0.9725 


TOTAL PRCNT 
. L E. 1 


1*193 

0.970* 

0.999* 

9889 

0.7515 

0.9913 

22 

0.3636 

0.8636 

5015 

0.6279 

0.97*3 

1258 

0.882* 

0.9968 

15682 

0.9333 

0.9999 

6031 

0.6201 

0.9833 

760 

0.7*7* 

1.0000 

19583 

0.7207 

0.99** 

**61 

0.6*99 

0.9662 

713 

0.7153 

0.9860 

7*098 

0.6950 

0.9*95 

151705 

0.7*78 

0.971* 







11. Figure Captions 


— L - 8‘ Target Area (TA) sampling using SAB method. Number 

of rejected TAs in each of the 40 latitude bands are 
shown by the dashed line. The solid line indicates the 
total number of TAs in the latitude band. Results are 
for the reflected flux. 


— B • — : Same as Fig. 1. The results are for the case of 

night-time emitted flux. 


— * S * — : a ) SAB, MLE 90° cutoff comparison. Differences 

of zonal 1 y averaged monthly mean albedo between the 
two methods. 

b) Linear regression relation between MLE with 
cutoff at 90° and SAB for albedo. The derived 
regression line parameters along with coefficient of 
correlation and the standard error of estimate are 
shown. Also shown is the value o where o 2 is 
the sum of the squares of the deviations. e 


Fig^ 4: a) Same as Fig. 3(a) for Day time LW flux 

differences . 

b) Same as Fig. 3(b) for Day time fluxes. 


Fi£i Is a) Same as Fig. 3(a) for Night time LW flux 
differences . 

b) Same as Fig. 3(b) for Night time fluxes. 


■ Fi R- a ) Same as Fig. 3(a) for Total LW flux 

differences. 

b) Same as Fig. 3(b) for Total fluxes. 


FigSj. Z^vlO: Same as Figs. 3-6 with the MLE method. 

Observations with satellite zenith angle greater 
than 75 are removed. 


— • ‘•8 s • — — Mercator maps of instantaneous albedo 

and Long wave flux fields (Day time, Night time and 
Total Long wave) between 70 S and 70 N. 

SAB method is used. 
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Figs. 15-18 : Mercator maps of instantaneous albedo and 

LW flux fields between 70 S and 70 N. MLE method, 
with cutoff at 75 in satellite zenith and with 
Nimbus-7 models, is used. 


Figs . 19-22 : Mercator projection of ERB parameter difference 

fields between SAB and MLE with cutoff in satellite 
zenith at 75 . 


Figs . 23-26 : Monthly zonally averaged clear and cloudy 

scene percentages with the use of Nimbus-7 models. 
Completely cloudy - Fig. 23; Clear - Fig. 24; 
Partly cloudy - Fig. 25; Mostly cloudy - Fig. 26. 


Fig . 27 : Instantaneous albedo with different zenith 

thresholds for Day 152. Results with bispectral 
thresholds are shown for comparison. 

Fig . 28 : LW fluxes for Day 152 as a function of zenith 

angle cutoff. 


Figs. 29-30 : Scene Id. Reliability Index for the 12 

scenes at the 2o and la levels for June 1, 1979. 
The terms NCLA etc. are defined in Frame 6. 
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Figure 14. 
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Figure 18. 
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